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Abstract 

Unified exact average capacity results for i-branch coherent diversity receivers including equal-gain combining 
(EGC) and maximal-ratio combining (MRC) are not known. This paper develops a novel generic framework for 
the capacity analysis of L-branch EGC/MRC over generalized fading channels. The framework is used to derive 
new results for the Gamma shadowed generalized Nakagami-m fading model which can be a suitable model for 
the fading environments encountered by high frequency (60 GHz and above) communications. The mathematical 
formalism is illustrated with some selected numerical and simulation results confirming the correctness of our 
newly proposed framework. 

Index Terms 

Average capacity, diversity, equal-gain combining (EGC), maximal-ratio combining (MRC), correlated channel 
fading. Gamma shadowed generalized Nakagami-m fading. 

I. Introduction 

Equal gain combining (EGC) is of practical interest in 60 GHz communications because its perfor- 
mance is comparable to that of maximal ratio combining (MRC) but it offers a greater simplicity of 
implementation (see [1] for an extended discussion on EGC and MRC performance difference). Due 
to high data-rate and coverage requirements of current, emerging and future high-frequency (60 GHz 
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or above) communication systems, the average capacity (AC) analysis of these two diversity combiners 
(i.e., EGC and MRC) becomes an important and fundamental issue from both theoretical and practical 
viewpoints. 

In literature, there are several papers dealing with the average symbol error probability (ASEP) analysis 
of the diversity receivers (see for example [1] and the references therein). Advances over the last decade 
on the symbol error performance analysis of EGC and MRC diversity receivers in fading channels has 
accentuated the importance of the moment generating functions (MGF) as a powerful tool for simplifying 
the analysis of diversity receivers. For example, the following identity has been widely used to simplify 
the symbol error performance analysis of EGC and MRC diversity receivers in fading channels. 



where erfc (■) is the complementary error function [2, Sec. (6. 13)], and where 'jend is the total signal- 
to-noise ratio (SNR) at the diversity receiver. On the other hand, and to the best of our knowledge, 
published papers dealing with the AC analysis of EGC and MRC diversity combiners over fading channels 
have been scarce when compared to those concerning the ASEP performance [1]. In particular, Bhaskar 
derived in [3] the average capacity of L-branch EGC relying on the Gamma approximation of the sum of 
mutually independent and identically distributed Rayleigh random variables (RVs). In addition, using an 
MGF-based approach, Hamdi obtained in [4] a new expression for the average capacity of MRC diversity 
combiner over arbitrarily correlated Rician fading channels. More recently, Di Renzo et. al proposed a new 
framework in [5] in order to compute the average capacity of MRC diversity combiner over generalized 
fading channels through the medium of the exponential integral Ej transform. However, the MGF-based 
approaches developed in [4], [5] were limited to the capacity of calculation of MRC diversity receivers 
and are not easily extendible to the computation of the capacity of EGC diversity receivers. In this paper, 
we show that it is actually possible to express the conditional capacity log2 (1 + 'Jend) in a form similar 
to (1), which facilitates the development of a new unified MGF-based approach for the calculation of 
the ergodic capacity in arbitrarily correlated/uncorrelated fading channels. More specifically, we present a 
unified MGF based average capacity computation not only for the L-branch MRC diversity receiver but 
also for the L-branch EGC diversity receiver over a wide variety of fading channels and for an arbitrary 
number of diversity branches. 

The remainder of this paper is organized as follows. In Section 11, a unified capacity analysis of diversity 




(1) 
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receivers over generalized fading channels is introduced and some key results are presented. In Section III, 
after the introduction of Gamma-shadowed generalized Nakagami-m (GNM) fading channel model, the 
exact average capacities for the EGC and MRC diversity receivers over Gamma- shadowed GNM fading 
channels are derived and many special cases are deduced. Numerical examples are then given in Section IV 
to illustrate the mathematical formalism. Finally, the main results are summarized and some conclusions 
are drawn in the last section. 

II. An MGF-based Capacity Analysis of Diversity Combiners 

For EGC and MRC diversity receivers, before the signals on the diversity branches are being summed 
to form the resultant output, the signals on the diversity branches are first co-phased and then weighted 
equally in EGC or weighted with the fading envelopes in MRC. The instantaneous SNR 7e„d at the output 
of the diversity receiver can be generically written as 



"^end 



L \ 1 



where the parameters p E {1,2} and q E {1,2} are chosen as 

f (1,2), EGC 

{p,q) = { . (3) 

[ (2,1), MRC 

In (2), L denotes the number of branches and Es/Nq is the transmitted SNR per symbol, and for i E 
{1, 2, 3, ... , L}, TZg is the ith branch fading. 

Considering the (instantaneous) Shannon capacity of the diversity receiver (i.e., EGC or MRC) with 
bandwidth W over fading channels (i.e., C^^^^ = W^log2(l +7end)), the average ergodic channel capacity 
defined as Cavg = 'E>[W logg (1 + 7end)], where E [■] denotes the expectation operator, can be obtained by 
averaging the instantaneous capacity C^^^^ over the probability density function (PDF) of ^end, namely 

/•oo 

Cavg = W l0g2 (1 + 7) P7eW (?) d^, (4) 

where p^^^^ (7) is the PDF of the instantaneous SNR jend (that is, 'jend is generically defined in (2)). 
Due to several reasons (e.g., insufficient antenna spacing or coupling among radio frequency (RF) layers), 
correlation may exist among diversity branches of the receiver. With or without that, the average capacity 
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using (4) involves an L-fold integral given by 

oo oo oo , 

Cavg = W f f... ^log2 [l + ^^-^==( J]rf j jp7e^,7e^,...,7e^ (ri,r2,...,r^)rfridr2. ..rfr^, (5) 




where ^'7^J,7^2,...,7^^ ('"d ^"2' • • • ' ''^l) '^he joint multivariate PDF of TZ^, ■ ■ ■ , TZ^ fading envelopes. The 
L-fold integration in (5) is tedious and complicated in addition to the fact that it cannot be separated into 
a product of one dimensional integrals. In addition, it takes a long time to evaluate numerically, especially 
as the number of branches L increases. Thus, referring to (4), researchers in literature have tried to find 
the PDF of the instantaneous SNR ^end given in (2) in order to find the average capacity. Nevertheless, this 
technique is often complicated and tedious for generalized fading environment since it involves multiple 
convolutions / integrals even if the fading envelopes 7^^,7^2) • ■ • of the branches are assumed to be 
independent. Referring to (2), the Jensen's inequality [6, Eq. (12.41 1)], which is based on concavity of log 
function such that E [log2(l + 7end)] < log2(l+IE [7end]), and fractional moments [6, Eq. (1.51 1)], which is 
based on the infinite series of log2(l+7end) such that E [log2(l + 7end)] = - Y.n>i (-l)'^E [7e„J/ log(2"), 
are commonly used in particular to compute the AC approximately. The other mostly used way to compute 
the AC hinges upon the inverse Laplace transform (ILT) whereby the PDF of the instantaneous SNR 'jend 
can be approximated through the medium of applying the ILT on the MGF A^7^„^ (s) = E [exp {—s^end)]- 
It is pertinent to say here again that the AC computation of diversity combiners (especially for the MGF 
of EGC since it is often more difficult than that of MRC due to the fact that it may not be possible to 
obtain the MGF of EGC receiver in general fading environments) becomes more difficult, problematic, 
and perplexing as the number of branches (i.e., L) increases. 

In what follows, we present a new exact and unified MGF-based approach that overcomes the difficulty 
mentioned above, and offers a generic single integral expression for the average capacity of EGC and 
MRC diversity combiners over generalized fading channels. 

Theorem 1 (Average Capacity of the Diversity Combiners over Correlated Not-Necessarily Identically 
Distributed Fading Channels). The exact average capacity of L-branch diversity combiner over mutually 
not-necessarily independent nor identically distributed fading channels with a bandwidth W is given by 



C ^ 



■'avg 



log (2) 



POO 

/ C,(.) 
Jo 



d 



ds, (6) 
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where p E {1,2} and q G {1,2} and are selected based on (3), and where the parameter $p ^ is defined 
as $p_g = ^ Es/NoL^P~i~'^y^, and where the auxiliary function Cg (s) is given by 



C, (s) = -Rlj 



1 

si 



(1,1), (0,1) 



(7) 



where H™^"[-] represents the Fox's H function [7, Eq. (1.1.1)]^'^. Moreover, M^p (s) = E[exp(— s Yln'^'i)] 
is the joint MGF for the p-exponent oflZ = {Tii, 7^2, ■ ■ ■ , ^l) feeding envelopes of the branches. 

Proof: See Appendix A. ■ 
Note that, in order to find the average capacity, the proposed MGF-based technique in Theorem 1 
eliminates the necessity of finding the PDF of the instantaneous SNR 7end through the ILT of the joint 
p-exponent MGF M.^p (s). Shortly, Theorem 1 suggests that one can readily obtain the average capacity 
of the diversity receiver by using the joint p-exponent MGF M.^p (s). Additionally, the integral in (1) 
can be readily estimated accurately by employing the Gauss-Chebyshev quadrature (GCQ) formula [11, 
Eq. (25.4.39)], yielding 

r a 1 

(8) 



W ^ id 



which converges rapidly and steadily, requiring only few terms for an accurate result, where the coefficients 
Wn and Sn are defined as 

Sn = tan (5 cos (%^7r) + ^) and Wn = ^ / , (9) 

^ 2^ ^ 4iVcos2 (f cos(^7r) + f)' 

respectively, where the truncation index N could be chosen as = 50 to obtain a high level of accuracy. 
In addition, when there is no correlation between the fading envelopes TZ = {7^]^,7^2, • • • the 
branches of the diversity receiver, the average capacity is given in the following corollary. 

Corollary 1 (Average Capacity of the Diversity Receiver over Mutually Independent Not Necessarily 
Identically Distributed Fading Channels). The exact average capacity of L-branch diversity receiver over 
mutually independent and non-identically distributed fading channels with the bandwidth W is given by 



fc=i 



'For more information about the Fox's H function, tfie readers are referred to [7], [8] 

^Note that the Fox's H function is still not available in standard mathematical software packages such as Mathematical® and MaplexM. 
However, using [9, Eq. (8.3.2/22)], an efficient mathematica implementation of this function is available in [10, Appendix A]. 
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where, for i E {1,2,..., L}, J^-jiP (s) = E [exp {—sTZ^)] is the MGF of the fading 71^ that the ith branch 
is subjected to. 

Proof: When there is no correlation between the fading envelopes TZ = {71^, 7l2, ■ ■ ■ , 7?.^}, one can 
readily write M.j^p{s) = Y[e=i-^'R-'' (■^)' whose derivative with respect to s is given by 



d_ 

ds 



M-uJi {s) 



(11) 



fe=i 



Finally, substituting (11) into (6) results in (10), which proves Corollary 1. ■ 
Despite the fact that the novel technique represented by Theorem 1 and Corollary 1 are easy to use, the 
numerical computation of the auxiliary function Cg (s) can also be done using the more familiar Meijer's 
G function, which is available in standard mathematical software packages such as Mathematica® and 
MaplcTM, as shown in the following corollary. 

Corollary 2 (Meijer's G Representation of the Auxiliary Function (s)). The auxiliary function Cg (s) 
can be given in terms of the more familiar Meijer 's G function as follows 

1,0 



-1 



-.G 



q{2 



TX 



1,2 

g+2,2 



2? 



1' 



(12) 



where h!"^! = -, 

(n) n ' n ' 



with X G C and ri G N. 



Proof: See Appendix B. ■ 
Let us consider the special cases (g = 1 for MRC and g = 2 for EGG) of the auxiliary function C^ (s) 
in order to check analytical simplicity and accuracy: 



Special Case 1 (Maximal Ratio Combining). For L-branch MRC diversity receiver (i.e., q 
auxiliary function Cmrc (s) = Cq {s)\^^-^ can be obtained as 



C 



MRC 



is) = -Gli 



'1 


1,1" 


s 






1), the 



(13) 



by means of applying [9, Eq. (8.2.2/9)] on (12). Utilizing [9, Eq. (8.2.2/14)] and [9, Eq. (8.4.11/1)] 
together, (13) reduces further to 

CMRcis)=Ei{-s), (14) 
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where Ei(-) is the exponential integral function [2, Eq. (6.15.2)]. ^. Then, referring to Theorem 1, the 
average capacity of the L-branch MRC receiver can be given by 

W 



nMRC _ 

log (2) 



/•oo 






/ Ei(-.) 






/o 


\N, J. 



ds 



(15) 



which is in perfect agreement with [5, Eq. (7)] when choosing the transmitted power is unit (i.e., 
Es/Nq = 1). In addition, when the branches are subjected to mutually independent and non-identical 
fading distributions, the average capacity C^y^^ can be also given, referring to Corollary 1, as follows 



C 



MRC 
avg 



w 



log (2) 



L r 



Eil 



1=1 



ds \No 



k = l 



(16) 



Special Case 2 (Equal Gain Combining). Note that, referring (7) with g = 2, the auxiliary function for 
L-branch EGC diversity receiver, i.e., Qegc {s) = Cg {s)\^^^ can be re-written as 

by means of setting q = 2 in (12). Eventually, using [9, Eq. (8.4.12/4)], the auxiliary function for L-branch 
EGC diversity receiver Cegc {s) simplifies to 



Cegc (s) = 2Ci {u) 



(18) 



where Ci (x) is the cosine integral function [11, Eq. (5.2.27)]^. Then, using Theorem 1 with (18), the 
average capacity of the L-branch EGC diversity receiver can be readily expressed as 



C. 



EGC 
avg 



2W 



Ci (u) 



9 ^A 



Es , 

iVnL' 



ds 



(19) 



in general when the fading T^i, 7^2, • • • , are subjected to are correlated. When the branches are 



subjected to mutually independent and non-identical fading, the average capacity Cj^^*" can be also given. 



referring to Corollary 1, as 

r,EGC ^ 

log (2) 



/"OO 

/ Ci(«)^ 



^ KA 



Es 

iVnL' 



A,— 1 



avg 



NnL 



■s ds. 



(20) 



In the following section, the model of Gamma shadowed GNM fading channel will be introduced and 

^Note that both the cosine integral function Ci (x) = — cos (t)/tdt and the exponential integral Ei (x) = — exp {—t)/tdt are 
implemented as a built-in function in the more popular mathematical software packages such as Mathematica® and MaplexM. 



IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. X, NO. XX, DEC. 2010 



8 



then the average capacity of both L-branch MRC and L-branch EGC diversity receivers will be derived 
for Gamma- shadowed GNM fading channels. 

III. Average Capacity of Diversity Combiners over Gamma-Shadowed Generalized 



As an example for the application of both Theorem 1 and Corollary 1, we assume that Gamma 
shadowing affects the GNM [10] fading channels, that is, the local mean power of the fading is a Gamma 
RV. For example, in 60 GHz non-line-of-sight (NLOS) propagation, standard deviation of shadowing is 
typically larger than that of propagation at 5 GHz. It is therefore not a misstep to assume that the local 
mean power of channel fading is a RV distributed over (0, oo). Hence, the Gamma- shadowed GNM fading 
model can be accommodated to pretest the evaluation of different wireless communications in 60 GHz 
non-line-of-sight (NLOS) environment. Thus, we first derive the analytical expressions for both PDF and 
MGF of the fading amplitudes a = {a-^^, a^-, • • • , otj^ for the branches of L-branch EGC in a Gamma- 
shadowed GNM fading channel. Using these results, we will find the exact average capacity for the EGC 
over a Gamma- shadowed GNM fading channel, and will give as examples of the simplified expressions 
for different special cases. 

We first derive the analytical expressions for both PDF and MGF of the fading envelopes 7? = 
{7^]^, 7^25 ■ ■ ■ '^he branches of L-branch diversity receiver in a Gamma-shadowed GNM fading 

channel. Using these results, we will find the exact average capacity for the diversity receiver over a 
Gamma-shadowed GNM fading channel and enumerate the different special cases. 

A. Gamma-Shadowed Generalized Nakagami-m Fading Channels 

Let us consider L > 1 mutually independent and non-identical GNM RVs each representing 

the fading amplitude the L-branch diversity combiner is subjected to and each having the PDF 



where the parameters > 1/2, > and fi^ > are the fading figure, the shaping parameter and 
the local mean power of the ith GNM RV and = F (m^ + 1/^^) /F (m^). Furthermore, the special or 
limiting cases of the GNM distribution are well-known in literature as the Rayleigh (m^ = 1,^^ = 1), 
exponential (m^ = 1,,^^ = 1/2), Half-Normal (m^ = 1/2,,^^ = 1), Nakagami-m (,^^ = 1), Gamma 
(^^ = 1/2), WeibuU (m^ = 1), lognormal (m^ oo,^^^ 0), and AWGN (m^ oo,,^^ = 1). 



NAKAGAMi-m Fading Channels 




(21) 
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As mentioned before, let the local mean power, of the GNM fading amplitude for the ith branch of 
the L-branch diversity combiner has, due to the shadowing, distribution with the Gamma PDF given by 

where fi^,^ is the average power of shadowing in the area of interest, and where m^^ inversely reflect the 
shadowing severity such that the severity of the shadowing decreases as the value of m^^ increases. For 
example, in the limit case m^^ — oo, the distribution of the local mean closes to the Dirac's distribution 
as pq^ = 5{^l — i^gi). Hence, there is no shadowing effects. 

Eventually, averaging (21) with respect to fi^, i.e., Pa^ {ci)Pn^ (^^) dVLf^, then utilizing [7, Theorem 
2.9] with [7, Eq. (2.9.4)], we obtain the PDF of the Gamma- shadowed GNM fading as introduced in the 
following definition. 

Definition 1 (Gamma-Shadowed Generalized Nakagami-m RV). The distribution TZi follows an Gamma- 
shadowed GNM distribution if the PDF of Tie is given by 



r(m,jr(m^) V / V it ^si iJ 

where (0.5 < < oo) and (0 < .^^ < oo) represent the fading figure (diversity severity / order) 
and the shaping factor, respectively, while rrigf^ (0.5 < m^^ < oo) and fl^^ (0 < < oo) denote the 
severity and the average power of shadowing, respectively. Moreover, F (■,■,■, •) is the extended incomplete 
Gamma function defined as V (a, x, 6, (3) = r""^ exp (— r — br~^) dr [12, Eq. (6.2)], where a and x 
are complex parameters, f3 > and b is a complex variable. 

In what follows, the shorthand notation TZ ~ J\fs{m,^,mg,^lg) denotes that TZ follows a Gamma- 
shadowed GNM RV with the fading figure m, the shaping parameter ^, the shadowing severity and 
the shadowing average power 

Let us consider some special cases of (23) in order to check validity. In fact, this PDF is a very 
general shadowed PDF which includes many special cases as explained in the second paragraph of this 
section. For example. By using [11, Eq. (6.1.47)] with the Mellin-Barnes contour integral representation 
[7, Eq. (1.1.1)] of (23), the PDF reduces to the PDF of the GNM distribution [10, Eqs. (1) and (2)] 
for rrigg — )■ oo as expected. Here again, by using [12, Eq. (6.41)], the PDF is reduced into the PDF of 
Gamma-shadowed Nakagami-m distribution [13, Eq. (9)] when the shaping parameter ^£ = 1. Furthermore, 



IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. X, NO. XX, DEC. 2010 



10 



the special or limiting cases of the Gamma- shadowed GNM distribution are well-known in literature as 
exponential- shadowed Rayleigh (m^^ = l,m^ = 1,,^^ = 1), K distribution (m^^ = 1,,^^ = 1), generalized- 
K distribution {^^ = 1), Rayleigh (m^ = 1,^^ = 1, m^^ -> oo), exponential (m^ = 1,^^ = 1/2, m^^ — )■ oo), 
Half-Normal (m^ = 1/2,^^ = l,m^^ — > oo), Nakagami-m (^^ = 1,772^,^ — )■ oo). Gamma (^^ = 1/2, m^^ — )■ 
oo), WeibuU (m^ = l,m^g -> oo), lognormal (m^ -> oo,^^ — )• 0,m^^ — )■ oo), and AWGN (m^ — > 
oo,^^ = l,m^^ — oo). Regarding these special and limit cases, the Gamma-shadowed GNM distribution 
TZ ~ J\fs{mi, ^ 1,711^^,^1 ^(^) has the advantage of modeling the envelope statistics of most known wireless 
and optical communication channels. Accordingly, it provides a unified theory as to model the envelope 
statistics. 

Referring to Theorem 1 and Corollary 1, we need to obtain the MGF of the fading envelope TZe ^ 
J\fs{m^,^^,m^^,il^^), i.e., M-jZc i-s) = E [exp {—sTZe)] for 3ft {s} G M"*" in order to find the exact average 
capacity of EGC, and also we need to obtain the MGF of the fading power 7^ = TZj, i.e., Ai-jij (s) = 
E [exp {—sTZj)] for 3ft {s} G in order to find the exact average capacity of MRC. As such, in the 
following theorem (i.e.. Theorem 2), these MGF functions are obtained in a unified closed form such that 
we can readily reduce it to the MGF of the £th branch of EGC and for that of MRC when the values 
p = 1 and j> = 2 are selected, respectively. 

Theorem 2 (Unified MGF of Gamma-Shadowed GNM RV). The unified MGF of the Gamma-shadowed 
GNM envelope distribution IZ^ ~ J\fs{m^, m^^, i.e., M-np (s) = E [exp {—slZ\)] is given by 



M-RV {s) 



r(m,^)r(m^; 



tt2,1 
■J^l,2 



n 



si 



(1,2) 

{m^^,p),{m^,p/i^ 



(24) 



with the convergence region 3ft {s} G 



Proof: Note that the unified MGF M-jiv (s) = E [exp {-sTZf}], Tie ~ J^sirn-e, it^ '"^s^, ^st) ^^''^ readily 
be given as A^t^p (s) = exp (— sr^) Pt^^ (r), where substituting the Fox's H representation of both 
extended incomplete Gamma function (i.e., [12, Eq. (6.22)]) and exponential function [7, Eq. (2.9.4)] 
results in 



r(m,^)r(mj j r 



tt1,0 

: 1^0,1 



sr'' 



(0,1) 



H 



2,0 
0,2 



{m^^,l),{m^,l/i^ 



dr. 



(25) 



Eventually, applying [7, Theorem 2.3] on (25), one can readily obtain the MGF of the Gamma- shadowed 
GNM distribution given in (24), which proves Theorem 2. ■ 
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Now, let US consider some special cases in order to check analytical simplicity and accuracy of (24). 
When setting the shadowing severity m « — )■ oo, and applying lim ^ 1, where 1 61 <^ a and 

|c| <C a, on the Mellin-Barnes integral representation [7, Eq. (1.1.1)] of (24), the unified MGF is, as 
expected, reduced to the MGFs of GNM [10, Eq. (2)] and generalized Gamma [14, Eq. (11)] for the 
values p = 1 and p = 2, respectively. 

Note that the unified MGF given by (24) may lead to some computation difficulty to compute due to the 
fact that the implementation of the Fox's H function is currently not available in standard mathematical 
software packages but an Mathematica® implementation of this function is offered by the authors in [10, 
Appendix A]. As such, it may be useful to represent (24) in terms of Meijer's G function with the aid of 
[9, Eq. (8.3.2/22)]. More specifically, (26) is the Meijer's G representation of (24) for the rational values 
of the parameter ^ (that is, we restrict ^ to ^ = k/l, where k and / are arbitrary positive integers.), namely 



. . f , y^{kp)"'^^{lp)"'' ^ r^kp+ip,2k 

J^^-RF \S) — kp.ip,, „ ^2k,kp+lp 

(2vr)^+^+'^-^rK)r(m,,) 



(?) 



kp / 1 \lp 



(kpnip) 



_^{-2k) 
-{2k) 
„(m^^) „(m^) 
-{kp) ^-{Ip) 



(26) 



It may be useful to notice that the rational representation of E M"*" remains essentially unchanged if 
there does not exist any rational number close enough to fulfilling the condition \k/l — < e//^, 
with e chosen to be 10^^. For more accuracy to rationalize the conditional parameter e can be chosen 
much smaller. Nevertheless, the number of coefficients of the Meijer's G function in (26) gets higher as e 
gets smaller, so much so that its computation efficiency considerable reduces and its computation latency'* 
increases. In consequence, the Fox's H function is preferable in this case since its computation efficiency 
is much better than that of Meijer's G function in this situation. 



B. Unified Average Capacity of Diversity Combiners 

Let us find the derivative of the unified MGF given by either (24) or (26) with respect to s since we 
need it in order to find the average capacity of L-branch diversity combiners operating over Gamma- 
shadowed GNM fading channels. Referring to the relation with an MGF Ai-jz^ (s) and its derivative, i.e., 
§l-M.ni {s) = — E [7?.^exp {—sTZ^)], the derivation of the unified MGF for diversity combiners operating 
over Gamma- shadowed GNM fading channels in give in the following theorem. 

Theorem 3 (Derivative of the Unified MGF of Gamma-Shadowed GNM RV). The derivative of the unified 

''The computation latency of Meijer's G function G^," [■] is primarily addressed by the total number of coefficients p + q. 
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MGF for the Gamma-shadowed GNM envelope distribution TZg ~ N's{rn^, m^^, i.e., ^Ai-jiP (s) 
— E [T^^exp {—sTZ^)] is given by 



ds ' 



(1,2), (0,1) 
[1,1), {m^^,p), {m^,p/Q _ 



(27) 



with the convergence region 3? {s} G 



Proof: Using either [7, Eq. (2.2.2)] or [9, Eq. (8.3.2/15)], the proof is obvious. ■ 
Again, following the same steps in the derivation of (26), (27) can be represented on the basis of the 
Meijer's G function for the rational values of the parameter ^ = k/l, where k and / are arbitrary positive 
integers. Accordingly, (27) can be given by 



d_ 

ds 



M-n?: {s) 



y^Al/^{kpr^^{lpr^-^k 



(27r 



pkp+lp+l,2k 
'^2k+l,kp+lp+l 



2fc \ 2fc / /3„m 



n 



kp 



{kpf'ilp)'" 



.{-2k) 
-{2k) ■■ 







1 ■^('^rf) 



(28) 



Finally, by employing Corollary 1 with (24) and (27), new exact single-integral expressions for the 
evaluation of the average capacity Cavg of L-branch diversity combiners over Gamma- shadowed GNM 
fading channels are immediately written as 



avg 



°° (-S) ^ j^3,l 
^ 1=1 



(f)2 o2 

^p,q^ 



(1,2), (0,1) 
[1,1), {m^^,p), {m^,p/ Q _ 



k = l 



(52 c2 



(1,2) 

(m,;-,p),(m;,,p/4; J 



ds, (29) 



where both the coefficient g and the auxiliary function Cg (s) are defined in Theorem 1. Furthermore, 



the coefficient Qi^ is defined as 



log(2) 



nLi r(m,^)r(m^ 



Additionally, referring to (8) (i.e.. 



by changing the variable of the integration in (29) as s = tan(6') and then using GCQ formula [11, 
Eq. (25.4.39)]), we specifically get a finite (A^-terms) sum approximation converging rapidly and steadily 
and requiring few terms for an accurate result as shown 



N 



^avg ~ 



Wr. 



Cq (Sn) 



n=0 



n 



si 
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(f)2 o2 
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l ^k^sk 
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(30) 



where Wn and s„ are defined in (9). In the sense of both that either special or limit cases of Gamma- 
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shadowed GNM A/'5(m^, m^^, model are commonly used fading models in the literature, and that 
the auxiliary function Cg (s) is a unified expression for diversity combiners (e.g., q = I for MRC and 
q = 2 for EGC), it is sufficient to show that the average capacity given by (29) is a unified expression not 
only for commonly used channel fading models but also for the commonly used MRC and EGC diversity 
combiners. For example, referring (14) (i.e., using the MRC special case of the auxiliary function Cg (s) 
given by (7)) and performing algebraic manipulations [7, Eqs. (2.1.1)-(2.1.5)] after choosing p = 2 for 
MRC, it is straight forward to show that, the unified average capacity (29) reduces to the average capacity 
of L-branch MRC diversity combiner over Gamma- shadowed GNM fading channels, namely 



-•MRC 

^avg 



)L+1 



Eil 



(1,1), (0,1) 
;i,l),(m^^,l),(m^,l/^^ 



X 



iVo/3,m,, 



(1,1) 

"^sfc,l),("^fc,l/4). 



ds, (31) 



Substituting the fading figures — > m, fading shaping factors = 1, the shadowing severities m^^ — )> oo 
and shadowing powers fi^^ = f2 for all ^ G {1, 2, . . . , L} in (29), and then using [7, Eqs. (2.1.1), (2.1.2), 
and (2.9.1)], the average capacity given by (31) reduces to the average capacity over mutually independent 
and identically distributed Nakagami-m fading channels, 

L-l 



c 



MRC 
avg 
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log(2)r^(m) 



Ei(-g) ^2,i 



2,2 



Nom 
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m 



ds. 



(32) 



Subsequently, note that we have G\1[u\ I] = m'' (1 + M)"T(a) [9, Eq. (8.4.2/5)] and GI'I 



u 



1,0 
l,a 



-u'' (1 + u) 



-a-1 



r(a) [9, Eq. (8.4.49/14)]. The average capacity of MRC diversity can be then attained 



in closed-form through the instrumentality of the Ei-transform equality Jg°°Ei(— n) (1 + an) ^ du 

-i/r(6)G 
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\ that is, 
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""''^ ~ log(2)r(mL + 1)^'2 



' 1_ 


0,0, 




m 


0, 


-1 



(33) 



where 7 = EsVl/Nq is the average SNR recovered by one branch of the MRC diversity combiner. Note 
that (33) represents an alternative compact closed-form expression (that is not limited to integer values of 
the fading figure m) to the result presented in either [15, Eqs. (19) and (20)] or [16, Eqs. (23) and (24)]. 

In addition, by choosing p = 1 and g = 2 for EGC, the auxiliary function Cg (s) simplifies into (18) 
and the unified average capacity given by (29) reduces to the average capacity of L-branch EGC diversity 
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combiner over Gamma-shadowed GNM fading channels, that is, 



nEGC 
avg 
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(1,1), (ii) 



ds, (34) 



which is obtained after using [7, Eqs. (2.1.1)-(2.1.5)]. Eventually, substituting the fading figures — )■ m, 
fading shaping factors = 1, the shadowing severities m^g = and shadowing powers fi^^ = Vl for all 
£ G {1, 2, . . . , L} in (34), and then using [7, Eqs. (2.1.1), (2.1.2), and (2.9.1)], the average capacity given by 
(34) reduces to the average capacity over mutually independent and identically distributed generalized-K 
fading channels as 

L-l 

pEGC 

^avg ~ I — r / „ "^2,2 



TT^ Jo 







p2,2 
^2,2 


ALNommg 


1^1 " 


E^n^s'^ 




E^n^s^ 





ds. 



(35) 



which can be readily computed by means of GCQ rule as seen in (8). 

It might be useful to notice that Tables I-III offer simplified expressions for the unified MGF and its 
derivative, for the variety of commonly used generalized fading channels in order to facilitate for the 
readers the use of our average capacity results for both MRC and EGG. 



IV. Numerical Results 

In this section, we provide some selected numerical results for the previous example, illustrating the 
average capacity of L-branch diversity receiver over Gamma- shadowed GNM fading channels. As seen in 
all figures (from Fig. 1 to Fig. 4), MRC gives better capacity/performance than EGG as expected, however 
it is a complex technique since it requires the envelope estimation of the channel fading. In addition, the 
minimum difference between their performances is obtained by two-branch combining. 

In Fig. 1, the average capacity of diversity receiver is depicted with respect to SNR (i.e., Es/Nq) 
for difference number of branches with Gamma shadowed GNM fading parameters V£ G {1,2,. 
rrig = 2,^^ = 2,m^^ = 3 and fi^^ = 1. This figure also displays the capacity per unit bandwidth (i.e., 
W = 1). Increasing the number of branches, i.e., L ^ 2, the average capacity increases but note that, 
regarding the relation among diversity gain and number of antennas, the diversity gain obtained by adding 
an antenna/branch decreases as the total number of antennas/branchs L increases. Note again that selected 
numerical and simulation results are in perfect agreement. 
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Amount of fading (AF) is another important statistical characteristic of fading channels, particularly 
in the context of applying diversity techniques in the transmission of the signals from transmitter to the 
receiver such as relay technologies. Shortly, this AF is associated with the fading figure / diversity order 
of the PDF given by (23), and for 7^^ ~ Afs{m^, m^^, fi^^), it is defined as = 
where Var [■] and E [•] are the variance and the expectation operators, respectively. As seen in Fig. 2, 
note that the large diversity gain is obtained by increasing fading figure / diversity order from 0.5 to 2.0. 
For example, a relay between transmitter and receiver or adding one more antenna to the transmitter may 
increase the diversity order from 1.0 to almost 2.0. For ^ 2, increasing the fading figure gradually and 
linearly increases the average capacity. In other words, increasing the number of relays or the number of 
the antennas at the transmitter more than 2 gradually and linearly increases the average capacity. Finally, 
note again that numerical and simulation results are in perfect agreement. 

When the signal recovered by the ith branch of L-branch diversity receiver from the wireless channel 
(i e {1,2, ...,L}) is composed of clusters of a multipath wave, each of which propagates in non- 
homogeneous environment^ such that they possess similar delay times and with the delay-time spreads 
of different clusters and their phases are independent [17]-[19]. In this case, the envelope of the received 
signal, i.e., 7^^ is considered as a non-linear function of multipath components. More specifically, let Xi 
and yi be the in-phase and quadrature Gaussian elements of the signal recovered from the Eth branch 
of the L-branch diversity receiver. Then, the envelope is represented as + ij^^l^^, where i = y/—T is 
the imaginer number and where is the shaping factor. Shortly, the shaping factor is sometimes not a 
sufficiently qualified parameter to comprehend and contemplate the fading conditions in some wireless 
communication applications. As such, for the PDF pn^ (r) of the fading envelope 7^^, the tail properties, 
i.e., both Umr^oo [pn^ ir)~\ and lim^^oo [dpn^ (r) /dr~\ axe changed by the shaping factor As seen in 
Fig. 3, the average capacity goes to zero when the shaping factor ^ goes to zero (V£ G {1,2,. ..,L}, 

= because, for < ,^ ^ 1, the tail properties approach to zero very fast with respect to the possible 
envelope values r G [0, oo), i.e., linv^oo [p7^^ (r)] = 0+ and Mnir^oo [dpiZf (r) /dr~\ = 0~. Also note that, 
for the higher values of shaping factor ^ ^ 1, the average capacity very gradually and linearly increases 
as seen in Fig. 3 as the shaping factor ^ increases. 

As mentioned at the beginning of the previous section, the link quality is affected by variation of 
the local mean power due to the shadowing caused by moving obstacles, scatters and reflectors between 

'Note that non-homogeneous wireless communications environment is very common in high frequencies such as 60 GHz or above due to 
the fact that the wave-length is very small when it is compared with the non-homogeneous (singular) environment. 



7^ 



/Var 



7^ 



IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. X, NO. XX, DEC. 2010 



16 



transmitter and the receiver. The intensity of shadowing on the branches of L-branch diversity receiver is 
characterized by the shadowing factors m^^ G [0.5, oo) for the branches £ G {1,2,...,L}. In Fig. 4, the 
average capacity is depicted with respect to shadowing factor (i.e., G {1,2,..., L}, m^^ = mj. 
Note that, as seen in Fig. 4, the average capacity does not change as the shadowing factor goes to 
infinity (i.e., — )• oo) since the variation of the local mean power diminishes as increases. 



In this paper, we presented a unified framework to compute the average capacity of diversity combining 
schemes (i.e., EGC and MRC) over fading channels. We also proposed a versatile fading model, which we 
term Gamma-shadowed GNM fading, in order to characterize the fading environment in high frequencies 
such as 60 GHz and above. Additionally, we derived novel closed-form expressions for the moment 
generating function (MGF) of both Gamma shadowed GNM fading and its special cases. Some selected 
simulations have been carried out for different scenarios of fading environment in order to verify the 
accuracy of the presented framework. Numerical and simulation results are in perfect agreement. 



This work was supported by King Abdullah University of Science and Technology (KAUST). 

Appendix A 
Proof for Theorem 1 

Note that, for q G {1,2} (i.e., g = 1 for MRC combining, and q = 2 for EGC combining), using the 
derivation equality Slog (1 + yX'') /dy = /{I + yR^), we can readily show that 



for n G M+. Using the equality z^'^ e^v{-sz) Ea^p{-yz'') dz = s""^/(s" + t/) [20, Eq. (5.2.3)], 
where E^^fs {■) is the Mittag-Leffler function [21, Eq. (1)], we get 



V. Conclusion 



Acknowledgments 




(A.1) 




(A.2) 
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Upon substituting ^ exp (—sX) = —X exp {—sX) into (A. 3), and then applying the well-known Leibnitz 
rule [11], it is easily shown that log (1 + X^) can be expressed as 



log(l + X« 



^exp(-sX) 



u 



-Eni \ — ]du 



u 



ds. 



(A.3) 



After substituting the Mellin-Bames representation of the Mittag-Leffler function, i.e., Ea^p {z) = <f^T (p) T (1 — p 
[21, Eq. (3)] and performing algebraic manipulations, (A.3) can immediately be expressed as 



log (1 + 



d_ 

ds 



exp (— sX) 



H 



1,2 
3,2 



1 



(l,l),(l,l),(l,g) 
(1,1), (0,1) 



ds, 



(A.4) 



by favor of the Mellin-Bames representation of Fox's H function [7, Eq. (1.1.1)]. Eventually, substituting 
(A.4) into (5) and using some algebraic manipulations, the average capacity of linear diversity receivers 
(EGC and MRC) can be readily given as in (6), which proves Theorem 1. 

Appendix B 
Proof for Corollary 2 

Note that, by means of [7, Eq. (1.1.1)], the auxiliary function Cq (s) given in (7) can be represented in 
terms of Mellin-Barnes integral as 



cjs)^-^^^-^±^m^$^s^^d. 



(B.l) 



27ri Jc r (1 - 2) r (1 + qz) 
with the convergence region 3ft {C} G ( — 1,0), where i is the imaginary number (i.e., i 



f^). Then, 

substituting Gauss' multiplication formula T (nz) = (27r)^^^""^ra'^^-^ HLi T + ^) [11, Eq. (6.1.20)] 
into (B.l) and using some algebraic manipulations, we get 



C,(3) 



-1 



r(i + z)r(-^)r(-^) fq 



si 



dz 



(B.2) 



Finally, applying [7, Eq. (2.9.1)] on the parenthesis part of (B.2), the auxiliary function Cg (s) can be 
derived as in (12), which proves Corollary 2. 
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Fig. 1. Average capacity versus tiie average power for different number of branches over Gamma-shadowed GNM fading channels 
Oil G {1, 2, . . . , L}, mg = 2,(^1 = 2, = 3 and fi^^ — 1). The number of samples for the simulation is chosen as N = 10000. 
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Fig. 2. Average capacity versus the channel fading figure for different number of branches over Gamma- shadowed GNM fading channels 
(Vi? G {1, 2, . . . , L}, 5^ = 2, rrigfi = 3 and il^^ = 1). The number of samples for the simulation is chosen as = 10000. 
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9.5 




Shaping Factor - ^ (V^ G {1, 2, . . . , L}, = 

Fig. 3. Average capacity versus the shaping factor for different number of branches over Gamma-shadowed GNM fading channels (V£ G 
{1, 2, . . . , L}, mg = 2, m^i = 3 and = 1). The number of samples for the simulation is chosen as A*" = 10000. 
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Fig. 4. Average capacity versus the shadowing factor for different number of branches over Gamma-shadowed GNM fading channels 
(it G {1, 2, . . . , L}, 7Ti£ = 2, = 2, rrigi = 3 and Q.^^ = 1). The number of samples for the simulation is chosen as A'' = 10000. 



TABLE I 

Unified MGFs of some Well-known Fading Channel Models 



Envelope Distribution, i.e., pn^ (r) 



Unified MGF M-j^p (s) and its derivative —M-jip{s), wliere tlie exponent p e {1,2} 



One-Sided Gaussian [1, Sec. 2.2.1.4] 
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defined over r 6 IR+, and where is the average power 
(i.e., > 0). Note that one-sided Gaussian fading coincides 
with the worst-case fading or equivalently, the largest amount 
of fading (AoF) for all Gaussian-based fading distributions. 
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where Gp^q [•] and Hp ^ [■] represent the Meijer's G function [9, Eq. (8.2.1/1)] and Fox's H function [9, Eq. (8.3.1/1)], respectively. In 

addition, the the coefficient hJ^! of the Meijer's G function is a set of coefficients such that it is defined as s!^! = — , . . . , 

with x g C and n 



Rayleigh [1, Eq. (2.6)] 
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defined over r S R+, and where is the average power 
(i.e., f2^ > 0). Note that Rayleigh fading distiibution has unit 
AoF (that is, AoF = 1). 
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Rayleigh distribution typically agrees very well with experimental data for mobile systems where no line-of-sight (LOS) path exists between 
the transmitter and receiver antennas [1, Sec. 2.2.1.1]. 



Nakagami-/M [1, Eq. (2.20)] 
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defined over r £ IR+, where is the average power, and 
where (0.5 < m^) denotes the fading figure. Moreover, 
r(-) is the Gamma function [6, Sec. 8.31]. 
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Note that the Nakagami-m distribution spans via the in parameter the widest range of amount of fading (AoF) among all the multipath 
distributions [1]. As such, Nakagami-^ (Hoyt) and Nakagami-n (Rice) can also be closely approximated by Nakagami-m distribution [I, 
Eq. (2.25)], [1, Eq. (2.26)]. 



WeibuU [1, Eq. (2.27)] 
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defined over r £ R+, where u)i = r(l + 1/^^) and where 
?f (0 < 5f ) denotes the fading shaping factor. Moreover, 
is the average power 
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where the Meijer's G representations are given for the rational value of the fading shaping factor (that is, we let = k/l, where k, 
and I ai'e arbitrary positive integers.) through the medium of algebraic manipulations utilizing [9, Eq. (8.3.2.22)]. In addition, note that if 
7?.^ is a sample of a Weibull distribution with the fading shaping factor then 7?." is also a sample of a Weibull distribution with the 
fading shaping factor ^a/o. 



TABLE II 

Unified MGFs of some Well-known Fading Channel Models 



Envelope Distribution, i.e., pn,, (r) 



Unified MGF M-j^p (s) and its derivative —M-j?p(s), wliere tlie exponent p e {1,2} 
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defined over r £ R+, where m^^ (0.5 < m^j.) is the fading 
figure, n^j, (0 < ^if.) is the average power, and 5^^, (0 < 
^^j.) is the accruing factor, of the fcth fading environment. 
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where r(-) is the Gamma function [11, Eq. (6.1.1)]. In addition. It may be useful to notice that the sum of the accruing probabilities 
k G {1, 2, . . . , of K possible fading environments is unit such that Ylk=i ^ik ~ ^■ 
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Nakagami-^ (Hoyt) [1, Eq. (2.10)] 
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defined over r 6 R+, where (0 < < 1) is the 
Nakagami-q fading parameter (that is, it is defined as ratio of 
the powers of the received signal's in-phase and quadrature 
with different standard deviations), and where f2j (0 < Q,p) 
is the average power. In addition, /q (•) is the zeroth order 
modified Bessel function of the first kind [11, Eq. (9.6.20)]. 
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where is defined as = 0.25(1 + q^)^ jq^, and is given by (q) = ^J^^^!;^ ((1 - + g^))^*", where fc e N. It may 

be useful to notice that the series expression of the unified MGF for the Nakagami-g (Hoyt) is converging very fast such that 10 summation 
teiTns is generally enough. 



Nakagami-« (Rice) [1, Eq. (2.15)] 
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defined over r £ R+, where (0 < n^) and (0 < C^) 
are the LOS figure and average power, respectively. 
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where .E^^. = r)^'' exp (— ) /fc! and the coefficient 9p = ^/2p/y'' (2Tr)^. In addition, the LOS figure i.e. rii is related to the Rician Ki 
factor by Ki = n| which con'esponds to the ratio of the power of the LOS (specular) component to the average power of the scattered 
component. 



K-Distribution [1, Eq. (2.15)] 
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defined over r 6 IR+, where m^^ (^ < m^^) denotes 
the shadowing severity, and (0 < f^^^) represents the 
average power Moreover, (■) is the nth order modified 
Bessel function of the second kind [11, Eq. (9.6.24)]. 
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It may be useful to notice that the shadowing effect in the channel disappears when m^^ approaches to infinity (m^ 
worst shadowing occurs when m^^ = ^. 
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00) such that the 



TABLE III 

Unified MGFs of some Well-known Fading Channel Models 
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Envelope Distribution, i.e., pm (r) 



Unified MGF Mfip (s) and its derivative — TMtjp (s), wliere tlie exponent p e {1,2} 
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Generalized-K [23, Eq. (5)] 
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defined over r 6 R+, where = m^^ + and = 
m^^—m^. Moreover, (0.5 < m^) and m^^ (0.5 < m^^) 
represent the fading figure (diversity severity / order) and the 
shadowing severity, respectively. (0 < f^^^) represents 
the average power. 
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where Qi = r(m^^)r(m^). It may be useful to notice that the shadowing effect in the channel disappears and generalized-K distribution 
turns into Nalcagami-m when m^^ approaches to infinity (m^^ — > oo) such that the worst shadowing occurs when m^^ = ^. 
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Composite Nakagami / Lognormal [1, Eq. (2.57)] M-j^p (s) 
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defined over r e M+, where ^^(m) = lot^"^? 
and where (dB) and cri(dB) are the mean and the standard 
deviation of channel shadowing. Moreover, (0.5 < m^) 
is the fading figure (diversity order), and (0 < Q^) 
represents the average power. 
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where, for n £ {1, 2, . . . , A^p}, {^^i„} and {a;„} are the weight factors and the zeros (abscissas) of the Afp-order Hermite polynomial 
[11, Table 25.10]. 



Composite Nakagami / Weibull [24, Eq. (4)] 
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defined over r £ R+, where (0 < H^) is the average 
power and (0 < 5^) denotes the Weibull (fading shaping) 
factor chosen to yield a best fit to measurement results. In 
addition, uj£ = T{1 + 1/$^) and m^^ (0.5 < ti^^) is the 
shadnMnng spverity 
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d. 

Note that Composite Nakagami / Lognormal distribution is the special case of Gamma-shadowed GNM distribution so the Meijer's G 
representation of the composite Nakagami / Lognormal distribution can be readily obtained by means of substituting = 1 and Q^^ = 
into both (26) and (28). 



Fox's H distribution [25, Eq. (3.1)], [26] 
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defined over r £ R+, and where /Cf S M and 
such two numbers that p-jz^ (r) dr = 1. 
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